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We study large-scale kinematic dynamo action due to turbulence in the presence of a 
linear shear flow, in the low conductivity limit. Our treatment is non perturbative in the 
shear strength and makes systematic use of both the shearing coordinate transformation 
and the Galilean invariance of the linear shear flow. The velocity fluctuations are assumed 
to have low magnetic Reynolds number (Rm) but could have arbitrary fluid Reynolds 
number. The equation for the magnetic fluctuations is expanded perturbativcly in the 
small quantity, Rm. Our principal results are as follows: (i) The magnetic fluctuations are 
determined to lowest order in Rm by explicit calculation of the resistive Green's function 
for the linear shear flow; (ii) The mean electromotive force is then calculated and an 
integro-differential equation is derived for the time evolution of the mean magnetic field. 
In this equation, velocity fluctuations contribute to two different kinds of terms, the "C" 
and "D" terms, in which first and second spatial derivatives of the mean magnetic field, 
respectively, appear inside the spacetime integrals; (iii) The contribution of the "D" terms 
is such that their contribution to the time evolution of the cross-shear components of the 
mean field do not depend on any other components excepting themselves. Therefore, to 
lowest order in Rm but to all orders in the shear strength, the "D" terms cannot give rise 
to a shear-current assisted dynamo effect; (iv) Casting the integro-differential equation 
in Fourier space, we show that the normal modes of the theory are a set of shearing 
waves, labelled by their sheared wavevectors; (v) The integral kernels are expressed in 
terms of the velocity spectrum tensor, which is the fundamental dynamical quantity that 
needs to be specified to complete the integro-differential equation description of the time 
evolution of the mean magnetic field; (vi) The "C" terms couple different components of 
the mean magnetic field, so they can, in principle, give rise to a shear-current type effect. 
We discuss the application to a slowly varying magnetic field, where it can be shown that 
forced non helical velocity dynamics at low fluid Reynolds number does not result in a 
shear-current assisted dynamo effect. 
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1. Introduction 

Large-scale magnetic fields in many astrophysical systems, from planets to clusters 
of galaxies, are thought to originate from dynamo action in the electrically conducting 
fluids in these objects. The standard paradigm involves amplification of seed magnetic 
fields, due to non mirror-s ymmetric (i.e. helical) turbulent flows, through the a-effect 
( Moffattl[l97i IParkeJlQTi ) . It is only relatively recently that the role of the mean shear 
in the turbulent flows is beginnng to be appreciated. Dynamo action due to shear and 
turbulence has received some a. ttention in the astro physical contexts of accretion disks 
( Vishniac fc Brandenbur j|l997t ) and galactic disks ( Blackman 1998f) . It has also been 
demonstrated that s hear, in coniunction with rotating turbulent convection, can dr ive a 
large-scale dynamo (jKapvla. Korpi fc Brandenberel 120081 : iHughes fc Proctojl2nn9h . We 
are interested in the more specific problem of large-scale dynamo action due to non- 
helical turbulence with mean s hear. Direct numeric al simulations now provide strong sup- 
port for such a shear dynamo. lYousef et al. ( 2008a ) demonstrated that forced small-scale 
non-helical turbulence in non-rotating linear shear fiows leads to exponential growth o f 
large-scale magnetic fields. These findings were later generalized bv lYousef et al. I (l2008bh 
to a shearing sheet mo del of a differentially rotat ing disk with a Keplerian rotation profile. 
The investigations of Brandenburg et all ( 20081 ) demonstrated the shear dynamo effect 
for a range of values of the Reynolds numbers and the shear parameter, and measured 
the tensorial magnetic diffusivity tensor. While the shear dynamo has been conclusively 
demonstrated to function, it is not yet clear what makes it work. This outstanding, 
unsolved problem is the focus of the present investigation. 

One possibility that has been suggested is dynamo action due to a "fluctuating a- 
effect" in turbulent flows which have zero mean helicities. In this proposal, large-scale 
dynamo action derives from the interaction of mean shear w i th fluctuations of helicity 



dVishniac fc Brandenbu^^ll997]: ISokolovlll997t IProctoil [2007i iBrandenburg et all [2008 



Rogachevskii fc Kleeorinll2008t ISchekochihin et al.ll2008 ). Another suggestion is that, if 



even transient growth makes non axisymmetric mean magnetic fields strong enough, 
they themselves rn ight drive motions which could lead to subcritical dynamo action 
( Rincon et al.ll2008|) . Yet anoth er possibility that has been suggested is the shear-current 
effect (jRogachevskii fc Kleeorin 2003, 2004, 2008) . In this mechanism, it is thought that 
the mean shear gives rise to anisotropic turbulence, which causes an extra component 
of the mean electromotive force (EMF), leading to the generation of the cross-shear 
component of the mean magnetic field from the component parallel to the shear flow. 
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However, there is no agreement yet whether the sign of suc h a couphng is favourable 
to the operation of a dynam o. Some analytic calculati ons (Radler fc Steoanov 20061 : 
Riidiger fc Kitchatinov 2006f) and numerical experiments ( Brandenburg et al. 2008 ) find 



that the sign of the shear-current term is unfavourable for dynamo action. A quasilinear 
theory of dynamo action in a linear shear f low of an incompressible fluid which has 
random velocity fluctuations was presented in Sridhar fc Subramanian ( 2009al lbl). Unlike 
earlier analytic work which treated shear as a small perturbation, this work did not 
place any restriction on the strength of the shear. They arrived at an integro-differential 
equation for the evolution of the mean magnetic field and argued that the shear-current 
assisted dynamo is essentially absent. It should be noted that the quasilinear theory 
of iSridhar fc Subramanian ( 2009a, b) assumes zero resistivity, and is valid in the limit of 
small velocity correlation times when the "first order smoothing approximation" (FOSA) 
holds. 

In this paper we present a kinematic theory of the shear dynamo that is non pertur- 
bative in the shear strength, but perturbative in the magnetic Reynolds number (Rm); 
this may be thought of as FOSA with finite resistivity. Thus we are not limited to the 
quasilinear limit of small velocity correlation times, and our conclusions are rigorously 
valid for velocity fiuctuations which have small Rm but arbitrary fluid Reynolds number. 
In § 2 we formulate the shear dynamo problem for small Rm. Using Reynolds averaging, 
we split the magnetic held into mean and fluctuating components. The equation for the 
fluctuations is expanded perturbatively in the small parameter, Rm. Using the shearing 
coordinate transformation, we make an explicit calculation of the resistive Green's func- 
tion for the linear shear flow. In § 3, the magnetic fluctuations and the mean electromotive 
force (EMF) are determined to lowest order in Rm. The transport coefficients are given in 
general form in terms of the two-point correlators of the velocity fluctuations. Galilean in- 
variance is a basic symmetry in the problem and is the focus of § 4. For Galilean invariant 
(G-invariant) velocity fluctuations, it is proved that the transport coefficients, although 
space-dependent, possess the property of translational invariance in sheared coordinate 
space. An explicit expression for the Galilean-invariant mean EMF is derived. We put 
together all the results in § 5 by deriving the integro-differential equation governing the 
time evolution of the mean magnetic field. Some important properties of this equation 
are discussed. In pa rticular, it is shown that , in the formal limit of zero resistivity, the 
quasilinear results of Sridhar fc Subramanian ( 2009al lbh are recovered. We also show that 
the natural setting for the integro-differential equation governing mean-field evolution 
is in sheared Fourier space. We prove a theorem on the form of the two-point velocity 
correlator in Fourier space; the velocity spectrum tensor and its general properties are 
discussed. We then express all the integral kernels in terms of the velocity spectrum ten- 
sor, which is the fundamental dynamical quantity that needs to be specified. Summary 
and conclusions are presented in § 6. 



2. The shear dynamo problem 

2.1. The small Rm limit 

Consider a Cartesian coordinate system with unit vectors (ei, 62, 63) erected on a comov- 
ing patch of a differentially rotating disk. Henceforth this will be referred to as the lab 
frame and we will use notation X = {Xi, X2, X^) for the position vector, and r for time. 
The fluid velocity is given by (5'Xie2 + v), where S is the rate of shear parameter and 
v{X,t) is a randomly fluctuating velocity fleld. The total magnetic field, S*°*(X,t), 
obeys the induction equation. 
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(-^ + SXi^ ] - S'Sf'ea = Vx (vxB'°') + rjV^B^°' (2.1) 

It is useful to note that the induction equation is unaffected by a uniform rotation of 
the frame of reference. So our coordinate system can refer to an inertial frame, or to 
a comoving patch of a differentially rotating disk. We study a kinematic problem in 
this paper, so will assume that the velocity field is prescribed. We also assume that the 
velocity fluctuations have zero mean {{v) = 0), with root-mean-squared amplitude t^rms 
on some typical spatial scale £. The magnetic Reynolds number may be defined as Rm = 
(frms^/??); note that Rm has been defined with respect to the fluctuation velocity field, 
not the background shear velocity field. To address the dynamo problem, we will use the 



approach of the theory of mean-field ele ctrodynamics ([Moffattl Il978t iKrause fc Radler 
19801: Brandenburg &: Subramanianll2005 ) . Here, the action of the velocity fluctuations on 
some seed magnetic field is assumed to produce a total magnetic field with a well-defined 
mean-field (B) and a fluctuating -field [b): 

_gtot ^ B + b, = B, {b) = (2.2) 

where ( ) denotes ensemble averaging in the sense of Reynolds. Applying Reynolds av- 
eraging to the induction equation (j2.ip , we obtain the following equations governing the 
dynamics of the mean and fluctuating magnetic fields: 

^ + SX,— )b - SBie2 = Vx5 + ijV^B (2.3) 

OT 0X2 I 



— + SXi-— ) b - Sbie2 = Vx (vxB) + Vx {vxb-{vxb)) + ?/V^6 



dr 'dX- 



(2.4) 

where £ = (vxb) is the mean electromotive force (EMF). The first step toward solving 
the problem is to solve equation (|2.4p for b, then calculate £ and obtain a closed equation 
for the mean-field, B{X, r). In the framework of the above mean-field theory, the shear 
dynamo problem may be posed as follows: under what conditions does the equation for 
B{X^t) admit growing solutions ? In particular, are growing solutions possible when 
the velocity field is non-helical (i.e. when the velocity field is mirror symmetric) ? 

The problem is, in general, a difficult one, but it can be approached perturbatively in 
the limit of small Rm. When Rm ^ 1, we can expand 6 in a series, 

b = + + + . . . (2.5) 

where b^"-* is of order b'"^^-* multiplied by the small quantity Rm. The equations gov- 
erning the time evolution of these quantitites are 

# + SX^^) - Sb["^e2 ^Vx{vxB) + TyV^fef") (2.6) 



dr 'dX- 



' d_ d 



SX,^\ 6(") - Sb^r^e2 = Vx (t;x6("-i) - (vxb^^-'^)) + jyV^b^ 
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for n ^ 1,2,... (2.7) 

Note that V X (i; X B) acts as a source term for b*-"-* , whereas the source term for b^"-* is 
Vx (vxb^"^"^^ - (v xb^'^'^^Yj . Once the b*^"^ have been determined, the mean EMF 
can be calculated directly by 

£ = {vxb) = (vx fb(°' + b(i) + 5(2) + _ ^\ (2.8) 



In this paper, we work to lowest order in Rni, so we need to work out only b^^'; equa- 
tion (12. 7p will not be used. 

2.2. The shearing coordinate transformation 

In this paper we will focus on the determination of the lowest order term. We also 
assume that the fluctuating velocity field is incompressible; i.e. V-t; = 0. Then the 
evolution of b'-^-' is governed by, 

^ + 5Xi-|^)b("' - Sb^°^e2 - {B-V)v - {vV)B + rjV^b'^''^ (2.9) 



dr 



We will now solve this equation for b'"-* and determine the mea n EMF. General method s 
of solving equations such as equation p.9p are presented in Krause fc Radleil ( IQSOl ). 
but we prefer to employ the shearing coordinate transformation because it is directly 
adapted to the problem at hand and greatly simplifies the task of writing down the 
Green's function solution. The {Xid/dX2) term makes equation (|2.9p inhomogeneous in 
the coordinate Xi. This term can be eliminated through a shearing transformation to 
new spacetime variables: 

Xi=Xi, X2=X2~ StXi , X3 = , t = T (2.10) 

Partial derivatives transform as 

d d ^ d d ^ St ^ ^ ^ ^ f2 11) 

dr dt ^ dx2 ' dXi dxi 8x2 ' 8X2 8x2 ' dX^, 8x3 

Therefore 



/ 9 8 
8Xp8Xp \ 8xp 8x2 



pp. pp. pp 

^ 2St-^^SH^^ (2.12) 



dxpdxp 8x18x2 8x2 

We also define new variables, which are component- wise equal to the old variables: 

H{x,t) = B{X,t), h{x,t) = b("^(X,T), u{x,t) = v{X,t) (2.13) 

Note that, just like the old variables, the new variables are expanded in the fixed Cartesian 
basis of the lab frame. For example, H — Hiei+H2e2+H3e3, where Hi{x, t) = Bi{X, r), 
and similarly for the other variables. In the new variables, equation (j2.9l) becomes, 



which can be expressed in component form as 

- J]V^] h,nix,t) = qm{x,t) (2.15) 



8t 
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where is given by equation ()2.12p . and 

qm{x,t) = [Hi - StSl2Hi]Uml - [ui - StdaUi] Hrnl + SSm2hi (2-16) 

We have used notation Umi = (dum/dxi) and Hmi — (dHm/dxi). Below we construct 
the Green's function for equation (|2.15p . 

2.3. The resistive Green's function for a linear shear flow 

Equation (|2.15p is hnear, homogeneous in x and inhomogeneous in t. Therefore, the 
general solution can be written in the form, 

hm{x,t) — j d^x' Grj{x — x' ,t, s) hmix' , s) 

^ Is ^^'/ ^^^'"^''(^ ^ ' foranys<t, (2.17) 

where Gn{x, t, t') is the resistive Green's function for the linear shear flow, which satisfies. 



^^-TjV^^G„{x,t,t') = (2.18a) 

lim G„{x,t,t') = S^{x) (2.18&) 

Grf{x,t,t') is non-zero only when < t' < t. (2.18c) 

G,j{x — x' ,t,tQ) = J d^x Gj^{x — x ,t,s)G,j{x — a;',s,io); for tg < s < t . (2.18rf) 
Let us define the spatial Fourier transform of the Green's function as, 

Gr,{k,t,t') = J d^x exp{-ik-x)Gr^ix,t,t') (2.19) 

where k, being conjugate to the sheared coordinate vector x, can be regarded as a sheared 
wavevector. Then 



dGr, 

dt 



TjK^{k,t)Gr, = 



(2.20) 



where, in equation (|2.20p . K^{k,t) — (fci — Stk2f' + A:| + /cf. It is now straightforward 
to write down the solution: 



Grj{k, t, t') — exp 



77 / dsK^{k,s) 



= exp 



-rj (^k^t - t') -Ski k2{t^ - t'^) + y fc2(^3 _ ^/3) 



(2.21) 



where, as per equation (|2.18cP above, t > t' . Note also that Grj{k,t,t') is a positive 
quantity which takes values between and 1, and that it is an even function of fc and k^. 

We now take the inverse Fourier transform of equation ()2.2ip to get Gri{x,t,t'). It is 
convenient to write this as 



Gri{x,t,t') — 



J (^ '^''P [ik-x - 7]{t - <'){fc2 + k, kj}] (2.22) 
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where Tij is a 2 x 2 symmetric matrix whose elements are given by, 
Tu — 0, Ti2 = T21 = 



(2.23) 



2 ' ^ " 3 

The integral in equation (j2.22l) can be evaluated by diagonalising the matrix Tij . It proves 
useful to express Grf{x,t,t') in terms of the principal-axes coordiates, x = {xi, X2, X3). 
These are defined by the orthogonal transformation, 



v^3 ; 



V 







cos( 







\ 



1 / 



X2 



V ^3 / 



(2.24) 



which is a time-dependent rotation of the coordinate axes in the X1-X2 plane. The angle 
of rotation, 0, is determined by 



tan^ = / + ^l + p 

1 S{e+tt' + t'^) 



f 



it + t') 



(2.25) 



Note that depends on the shear parameter, S, and the times, t and t' . Let us define 
the dimensionless quantitites. 



S 

1 - —(t + t') tan( 



1/2 



(T2 = 



1 + -(< + t') C0t( 



1/2 



Now we can write the Green's function as a sheared heat kernel, 



G^{x,t,t') = [AlTTjit ~ t')]~^^^ 



12 



r\2 



(t-f) 



X exp 



Ar]{t-t') 



7^:2 —2 



''2 I —2 



^1 



(2.26) 



1-1/2 



(2.27) 



which is equivalent to the one first derived in iKrause fc Radleii (jl97lh . 

We now note some properties of the Green's function. For convenience we choose the 
shear parameter, S, to be negative: then the quantities, / ^ 0, ^ 6* ^ 7r/2, cti ^ 1 
and ^ (72 ^ 1 . At fixed t and f , the Green's function is a Gaussian with long axis 
along xi, short axis along X2, and the intermediate axis along X3. To obtain some idea 
of the behaviour of the Green's function, it is useful to plot isocontours in the sheared 
coordinate space (xi, X2, cca) at different values of t and t'. Figure ^ displays isocontours 
in the X1-X2 plane at four different values of t for t' = 0; we have chosen X3 — and t' — 
in the interests of brevity of presentation. The figure is plotted in shearing coordinates, 
with respect to which diffusion is anisotropic and there is no advection. It may be noted 
that the Green's function shows a shearing motion against the direction of the actual 
shear. As t increases from zero to infinity, 9 (which is the angle the long axis makes with 
the a;i-axis) increases from 45° to 90°, and all the principal axes increase without bound. 




-40 -20 20 40 -40 -20 20 40 



Figure 1. Isocontours of the resistive Green's function Gr,{x,t,t') plotted in the xi-X2 plane 
of the shearing coordinate system, for t' = at four different values of t. Units are such that 
5* — —2 -jij — 1. Five isocontours at 90%, 70%, 50%, 30% and 10% of the maximum value are 
displayed. Panels (a), (b), (c) and (d) correspond to times t — 1, t = 5, t = 10 and t = 15. 



3. Magnetic fluctuations and mean EMF at small Rm 

3.1. Explicit solution for h{x,t) 

We are interested in the particular solution to equation (I2.15|) (i.e. the forced solution) 
which vanishes at t = 0. This can be written as 

K,{x,t) ^ J^dt' J d!^x'G,j{x~x',t,t')q„,{x',t') (3.1) 
Substituting the expression for from equation (|2.16p in equation p.ip . we have 
hm{x, = J d^x' Gjj{x - x', t, t') X 

X {[H'l - St'6i2H[] u'^i - [u'l - St'5nA] H'^i} 



Shear dynamo problem 



SSm2 dt' J d^x'G,,{x-x',t,t') hi{x\t' 
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(3.2) 



where primes denote evaluation at spacetime point {x'^t'). The solution is not yet in 
explicit form because the last term on the right side contains the unknown quantity 
hi{x',t'). Thus we need to work out the integral 

^ dt' J d^x' Gr,{x-x',t,t')hi{x',t') = dt' J d^x'G,,{x-x',t,t') X 

xj dt" J d^x" G,,{x' -x" ,t',t") X 



Hi - St SnH^ 



St 



where " means evaluation at spacetime point [x ,t ). Note that, on the right side, x' 
occurs only in the Green's functions. So, by using the property given in equation (|2.18(ij) . 
the integral over x' can be performed. Then 

/ dt' [ d^x'Gr,{x-x',t,t')hi{x',t')^ [ dt' [ dt" [ d^x" Gr,{x-x",t,t")x 
Jo J Jo Jo J 



{ 



Hi ~ St 5i2H, 



u, — St 



S12U1 



The double-time integrals can be reduced to a single-time integrals because of the fol- 
lowing simple identity. For any function f{x,t), we have 

f dt' f dt" f d^x" f{x\t")= f dt" f d^x" f{x\t") f dt' 
Jo Jo J Jo J Jt" 



= dt (t-t ) d^x f{x ,t ) 



dt'{t~t') / d^x' f{x',t') 



where in the last equality we have merely replaced the dummy integration variables 
{x ,t ) by {x',t'). Then we have 



dt' J d^x' Gr,{x-x',t,t') hi{x'X)^ dt'{t-t') J 



)hi{x'X)^ J dt'{t-t') J d^x'G^{x-x ,t,t)x 

X {[H^ - St'SnH'^] u'^i - [u'l - St'Snu'^] H'^i) 
Therefore the forced solution to equation (|2.15p can finally be written in explicit form as 



h^(x, t)^ J dt' J d^x' G^{x - x', t, t') + S{t - t')5^2u'u] X 

X [H'l - St'SM 

- j^'^t' j d'^' Gr,{x - X', t, t') [ij;^; + S{t - t')5ra2H'^i] > 
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X [u'l ~ St'5i2u[] (3.3) 
This gives the magnetic fluctuation to lowest order in Rm. 

3.2. Explicit expression for the mean EMF 

To lowest order in Rm, the mean EMF is given by 5 = (vxb^^'^^ = {uxh), where 

equation p.3p for h should be substituted. The averaging, ( ), acts only on the velocity 
variables but not the mean field; i.e. {uuH) = {uu) H etc. After interchanging the 
dummy indices {l,m) in the last term, the mean EMF is given in component form as 

dt' [ d^x' Gr,{x - x',t,t') \aii{x,t,x',t') + S{t ~ t')l3a{x,t,x' ,t') 



X [H'l - St'5n H[] 

dt' j d^x' Gr,{x - x' ,t,t')[rii„d{x,t,x' ,t') - St'6„,2r]iii{x,t,x' ,t')] x 
x[i/L + S{t~t')dnH[J 

(3.4) 

Here, {a ^ l3 ,r]), are transport coefficients^ which are defined in terms of the uu velocity 
correlators by 

aii{x,t,x',t') = eijm.{Uj{x,t)Uml{x',t')) 
I3ii{x,t,x' ,t') = €ij2 {Uj{x,t)uii{x\t')) 

rjimi{x.,t,x' ,t') = eiji{u.j{x,t)um{x' ,t')) (3.5) 

It is also useful to consider velocity statistics in terms of vv velocity correlators, because 
this is referred to the lab frame. By definition, from feauation l2.13p . 

Um{x,t) = Vm{X{x,t),t) (3.6) 

where 

Xi ^ xi, X2 = X2 + Stxi , = X3 , T ^ t (3.7) 

is the inverse of the shearing transformation given in equation (j2.10p . Using 

d d d 

oxi dXi 8X2 
the velocity gradient Umi can be written as 

Uml = = ( + St^IIT^] Vm = Vml + St5iiV„i2 (3.9) 

dxi \ dXi 8X2 J 

where Vmi = {dvm/dXi). Then the transport coefficients are given in terms of the vv 
velocity correlators by 

au{x,t,x',t') = e^j^[{vjiX,t)vn,iiX',t')) + St' 611 {vj{X,t)vm2{X',t'))] 
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^u{x,t,x\t') = e^j2[{vj{X,t)viiiX',t')) + St'Sn {vj{X,t)vi2{X\t'))] 

riirnl{x,t,x',t') = eijl {vj{X,t)Vm{X',t')) (3.10) 

where X and X' are shorthand for 

X = (.Ti , .T2 + , X3) , X' = {x[ ,X2 + St'x'i ,x'^) (3-11) 

Equation p.4p . together with p. 51) or p.lOp . gives the mean EMF in general form. X 
can be thought of as the coordinates of the origin at time t of an observer comoving with 
the background shear flow, who was at x at time equal to zero. Similarly, X' can be 
thought of as the coordinates of the origin at time t' of an observer comoving with the 
background shear flow, who was at x' at time equal to zero. 



4. Galilean invariant velocity statistics 

4.1. Galilean invariance of the induction equation 

The induction equation ()2.1|) for the total magnetic field — and also equations p.3p and 
(|2.4p for the mean and fluctuating components — have a fundamental invariance property 
relating to measurements made by a speci al subset of all observers, called comoving 
observers in Sridhar fc Subramanian (l2009allbh . A comoving observer translates with the 
velocity of the background shear flow, and such an observer can be labelled by the 
coordinates, ^ = (^1,^2,^3)1 of her origin at time r = 0. At any time r, the origin is at 
position, 

Xc{i,r) = + (4.1) 

An event with spacetime coordinates {X^t) in the lab frame has spacetime coordinates 
(X,f) with respect to the comoving observer, given by 

X = X X,(^,t), f = r-To (4.2) 

where the arbitrary constant tq allows for translation in time as well. 



Let 



B'°\X, t) , B{X, t) , f ) , v{X, f) 



denote the total, the mean, the fluctu- 



ating magnetic fields and the fluct uating velocity field, respectively , as measured by the 
comoving observer. As explained in lSridhar &: SubramanianI ( 2009al lbf). these are all equal 
to the respective quantities measured in the lab frame: 



B'°\X, f) , t) , f) , v{X, f) 



iX,T),B{X,T),b{X,T),v{X,T)] 

(4.3) 

It is proved in ISridhar fc SubramanianI ('2009a'.'D) that equations ([2TT|) . ([231) and (|2^ 
are invariant under the simultaneous transformations given in equations (|4.2p and (|4.3p . 
This symmetry property is actually invariance under a subset of the full ten-parameter 
Galilean group, parametrized by the five quantities {^i,^2,£,3,tq, S); for brevity we will 
refer to this restricted symmetry as Galilean invariance, or simply GI. 

It is important to note that the lab and comoving frames need not constitute inertial 
coordinate systems. One of the main applications of our theory is to the shearing sheets 
which is a local description of a differentially rotating disk. In this case the velocity field 
will be affected by Coriolis forces. The only requirement is that the magnetic field satisfies 
the induction equation ()2.1|) . 



4.2. Galilean-invariant velocity correlators 

We now explore the consequences of requiring that the statistics of the velocity fiuc- 
tuations be Galilean-invariant. We consider the n-point velocity correlator measured 
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by the observer in the lab frame. Let this observer correlate vj-^ at spacetime location 
(i2i,Ti), with Vj^ at spacetime location (i?2,T2), and so on upto vj^ at spacetime lo- 
cation {Rn,Tn). Now consider a comoving observer, the position vector of whose origin 
is given by Xc{^,t) of equation (j4.ip . An identical experiment performed by this ob- 
server must yield the same results, the measurements now made at the spacetime points 
denoted by (/li + n), ; (i?2 + Xc(^, tj), rs) ; . . . ; (H„ + r„), t„). If the 

velocity statistics is GI, the ?i-point velocity correlator must satisfy the condition 

{Vj^{Rl,Ti) . . .Vj^{Rn,Tn)) = {Vj^{Ri+Xc{$,.Ti),Ti) . . . Wj^ (H„ + Xc , r„ ) , r„ ) ) 

(4.4) 

for aU (fii, . . . i2„ ; Ti, . . . t„ ; |). 

In the low Rm limit, we require only the two-point velocity correlators, for which 

{v,{R,t)v,{R!,t')) = {v,{R + X,{i,T),T)v,{R' + X,{i,T'),T')) (4.5) 

for all {R, R' ,t,t' ,^). We also need to work out the correlation between velocities and 
their gradients: 

{v,{R,T)v,iiR\r')) = jL(y^^R^r)v,iR',T')) 

= ^ HR + X^it t),t) V, {R' + X,(C, r'), r')) 

= {v,{R + X,(C, r), t) v,i {R' + X,{t r'), r')) (4.6) 

We want to choose {R, R' , r, r', ^) as functions of {x, x' , t, t') such that we can use equa- 
tions (|4.5p and (j4.6p to simplify the velocity correlators in equation (|3.10p . We note that 
equations p. lip and (14.11) give 

X = X,{x,t), X' = X,{x',t') (4.7) 

It is therefore natural to choose 

T = t, t' = t' (4.8) 
Thus the velocity correlators we require can now be written as 

{v,{X,t)vj{X\t')) = {v,{X,{x,t),t)vj{X,{x',t'),t')) 

{v,{X,t)vji{X'X)) = {v,{X,{x,t),t)v,i{X,{x\t'),t')) (4.9) 
Comparing equation (|4.9p with equations (|4.5p and ()4.6p , we see that if we choose 

R^Xc{x,t), R' = Xc{x',t') (4.10) 
then equation (|4.9p . together with equations (|4.ip . (|4.5p and ()4.6p implies that 
{v,{X,t)vj{X',t')) = {v,{R,t)vj{R',t')) 

= {vi{R + X,(|, t),t)v, {R' + X,(|, r'), r')) 
= {v,{Xc{x + t),t)v,{X,{x' + ^, t'),t')) 

Similarly 

{v,iX,t)v,i{X',t')) = {v,{X,{x + ^,t),t)vji{X,ix' + i,t'),t')) (4.11) 
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Figure 2. Galilean invariance of the two-point velocity correlator given in equation (|4.5(l . O 
labels the observer in the laboratory frame who correlates the velocity fluctuation at location 
R at time r with the velocity fluctuation at location R' at a later time r'. A and A' label a 
comoving observer the origin of whose cooordinate axes is at ^ at the initial time, and who 
makes an equivalent measurement at the times r and r'. 



Now it is natural to choose 



(4.12) 



Then 






t] ,t]vi{ Xr i Z ,t' ,t' 



(4.13) 



Similarly, 



(vi{X,T)Vjl{X',T')) ^ ivi [X 



X — X 



t ,t v.ilX 



X — X 



1 " If 



= Qijlix - x',t,t') 
We note that symmetry and incompressibility imply that 

R^j{r,t,t') = R.j,{-r,t',t) 

Q,jj{r,t,t') = 



(4.14) 



(4.15) 
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4.3. Galilean-invariant mean EMF 

The transport coefficients arc completely determined by the form of the velocity corre- 
lator. Using equations (|4.13p and (j4.14p in equations (I3.10|) and noting the fact that the 
velocity correlators defined above are functions only of {x — x'), t and t' , we can see that 
the GI transport coefficients, 

au{x,t,x',t') = e,jjn[Qj^i{x - x',t,t') + St' SiiQj,n2{x ~ x' ,t,t')] 

%i{x,t,x' ,t') = eij2[Qjii{x ~ x',t,t') + St' SiiQ ji2{x - x\t,t')] 

rji„ii{x,t,x' ,t') = eijiRjm{x - x' ,t,t') (4.16) 

The transport coefficients depend on x and x' only through the combination, {x — x'), 
which arises because of Galilean invariance. We can derive an expression for the G- 
invariant mean EMF by using equations (I4.16|) for the transport coefficients in equa- 
tion p.4p . We also change the integration variable in equation p.4[) to r = x — x' . The 
integrands can be simplified as follows: 

auix,t,x',t')[H'i^St'6i2H[] = e,jm[Qjmi + St'5nQj^2][H'i-St'5nH[] 

= eij„iQjmi{r,t,t')Hi{x - r,t') 

%i{x, t, x', t') [H[ - St'5i2H[] = e,,2 [QjU + St' Sn Q,i2] [H'l - St'S^H'^] 

= eij2Qjii{r,t,t')Hi{x - r,t') 

[Vinii - St'S,n2rjiii]H'i^ = eiji[Rj,n{r,t,t') - St'Sm2 Rji{r,t,t')] Hi„^{x - r,t') 

[%ni2- St'5m2%i2]H'^,n = £-^2 <5a [i?j>„ (r , i, t') - St'6^2Rji{r,t,t')]Hi^{x - r,t') 
Define 

Cjmi{r,t,t') = Qjrni{r,t,t') + Sit-t')Sm2QMr,t,t') 

Dj^{r,t,t') = Rjra{r,t,t') - St' 5ra2Rji{r,tX) (4.17) 
The mean EMF can now be written compactly as 

E,{x,t) = ^rjm dt' J d?r Gn{r,tX)Cjmi{r,t,t')Hi{x - r,t') 

dt' J dV G„(r, t, t') [e,ji + Sit - t')Sne,j2] x 

xDj^{r,t,t')HiUx-r,t') (4.18) 

5. Mean field induction equation 

5.1. Mean-field induction equation in sheared coordinate space 

Applying the shearing transformation given in equations (|2.10p and p. lip to the mean- 
field equation p.3p . we see that the mean-field, H{x, t), obeys 

B H 

- SSaHi = (Vx£), -t- r^V^H, (5.1) 
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where 

We note that the divergence condition on the mean magnetic field can be written as 

f) H 

V-H = ^ Hpp - StHi2 - (5.3) 

ClAp 

It may be verified that equation (|5.ip preserves the condition V- -Ff = . We now use 
equations (|4.18l) and (15.21) to evaluate Vx5. 

d£ ( d d \ 



J dt' J d^r G^{r, t, t')C„ni {r, t, t') [iZ/^ - St Sp.H'i^] 



dt' J d^rGrj{r,t,t')Dj„,{r,t,t') [e^pgCgji + S{t - t!)iiie.,pqeqji\ x 

where il[ = Hi{x — r,t'). Expanding tipqtqjm — {Sij Smp — Sim 5jp), the contribution 
from the C term is 

(Vx£)f = J^t' J d^rGr,{r,t,t')[C,pi{r,t,t')~Cpa{r,t,t')]x 

X [Hip - StSp.Hl,] (5.4) 

Evaluating the D term is a bit more involved. Again, we begin by expanding tipqtqji = 
(Sij Sip - Sii S-jp). Then we get 

(Vx5)f = ^ dt' J d''rG^{r,t,t')Dpm{r,t,t') x 

X {H'ipra - StSpiH-2m + ^ {t - t')6i2 [H[p^ ~ StSpiH^^m] } 

- j^'^t' j d'^ G^(r, i, t') D,m{r, t, t') [H^^^ - St'Hi^J (5.5) 

The second integral vanishes because the factor in [ ] multiplying Dim is zero: to see this, 
differentiate the divergence-free condition of equation (|5.3p with respect to Xm- We can 
now use equations ((^ and ((^ to write (Vx5) = (Vx^)*^ + (Vx£")^. Substituting 
this expression in equation (|5.1[) , we obtain a set of integro-difFerential equation governing 
the dynamics of the mean-field, H{x, t), valid for arbitrary values of the shear strength 
S: 



dt 



SSaHi = T^V^H, + dt' J d\G^{r,t,t') [C^mi{r,t,t') - Cmu{r,t,t')] 
x[Him{x-r,t') - St5miHi2{x-r,t')] 



16 S. Sridhar and Nishant K. Singh 

+ dt' J <i\G.,-i{r,t,t')Djm{r,t,t') x 

X [H,jm{x - r,t') ~ St5jiHam{x - r,t') + 

+ Sit - t')S,2 {Hijm{x - r,t') - StdjiHi2rn{x ~ r,i')}] 

(5.6) 

We note some important properties of the mean-field induction equation (j5.6l) : 

(a) The Djm{r^t,t') terms are such that (Vx£)j involves only Hi for i— \ and i = 3, 
whereas (Vx5)2 depends on both H2 and Hi. The implications for the original field, 
B(X,t), can be read off, because it is equal to H{x,t) component-wise (i.e Bi{X,T) = 
Hi{x,t)). Therefore, in the mean-field induction equation, the "D" terms are of such a 
form that: (i) the equations for Bi or B3 involve only Bi or B3, respectively; (ii) the 
equation for B2 involves both Bi and -B2. 

(&) Only the part of Cim;('', t, t') that is antisymmetric in the indices (i, m) contributes. 
We note that it is possible that the "C" terms can lead to a coupling of different com- 
ponents of the mean magnetic field. To investigate this, it is necessary to specify the 
statistics of the velocity fluctuations. 

(c) In the formal limit of zero resistivity, rj 0, the resistive Green's function, 
G{x,t,t') S{x). Then the mean-field induction equation simplifies to 

f)H 

- SdaHi = / di' [C„„i(0,i,O-C™z(0,i,O] X 

ot Jo 

X [Hi,Jx,t') - St5,nlHi2{x,t')] 



f dt' D,„,{0,t,t')[H,j„,{x,t') - StdjiH,2m{x,t') 

Jo 

+ Sit - t')S,2 {Hi,^{x,t') - St5jiHi2.n{x,t')}] 



(5.7) 



which is identical to that derived in Sridhar fc Subramanian (|2009allbh 



5.2. Mean-field induction equation in sheared Fourier space 

Equation ()5.6p governing the time evolution of the mean field may be simplified further 
by taking a spatial Fourier transformation. Let us define 

H{k,t) = j d^xH{x,t)eyij>{-ik-x) (5.8) 

and the quantities, 

hmi{k,t,t') = [ d^rGr,ir,t,t')CM{r,t,t')exp{~ik-r) 



(5.9) 



Both Iimi{k,t,t') and Jj„i{k,t,t') are to be regarded as given quantities, because they 
are known once the velocity correlators have been specified. Taking the spatial Fourier 



J„n{k,t,t')= / d^rG^{r,t,t')Dj„,{r,t,t')exp{-ik-r) 
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transform of equation ()5.6|) , we obtain, 

- SSaHi^'f]K^H, + iK^j dt' [i,^iik,t,t') - i^aik,t,t')\ Hi{k,t') 

-hnKj f dt' Jjrnik, t, t') \H,{k, t') + S{t - t')S,2Hi{k, t' 

(5.10) 

where K{k,t) = (fci - Stk2, fe, fcg) and K'^ = = {ki - S'ifcs)^ + fc| + fc^ . Once the 
initial data, i?(fc,0), has been specified, equations (jS.lOp can be integrated in time to 
determine H{k, t). Whereas these equations are not easy to solve, we note some of their 
important properties: 

(a) Only the part of limi {k, t, t') that is antisymmetric in the indices {i, m) contributes. 

(&) The time evolution of H{k, t) depends only on H{k, t') for ^ i' < <, not on the 
values of H at other values of k. Thus each k labels a normal mode whose amplitude 
and polarisation are given by H{k,t), the time evolution of which is independent of all 
the other normal modes. 

(c) When we have determined H{k,t), the magnetic field in the original variables, 
B[X,t), can be recovered by using the shearing transformation, equation (j2.10p . to 
write {x,t) in terms of the lab frame coordinates {X,t): 

/d^k ~ 
.^H{k,t)exp{ik-x) 

= J ^fl-(fc,r)exp[iK(fe,r).X] (5.11) 

where we have used K- X = k- x. Thus B(X,t) has been expressed as a superposition 
of the normal modes, each of which is a shearing wave, whose spatial structure is given 

by 

exp[iK(fc,T).X] = exp[i{{ki- Stk2)Xi + k2X2 + k3Xs}] (5.12) 

For non-axisymmetric waves, k2 ^ and, as time progresses, the shearing wave develops 
fine-structure along the Xi-direction with a time-dependent spatial frequency equal to 

(fcl - ^2). 

5.3. The integral kernels expressed in terms of the velocity spectrum tensor 

We have derived the integral equation satisfied by the mean magnetic field, to lowest 
order in Rm; in sheared coordinate space it is given by equations (15.61) . and in sheared 
Fourier space it is given by equations (|5.10p . One can proceed to look for solutions if the 
integral kernels are known. This means that either the pair [Cimi{r,t,t'),Djra{r,t,t')] 

or the pair Iimi{k, t, t') , Jj,n{k, t, t') needs to be specified. Here we show that all these 

integral kernels can be expressed in terms of a single entity, which is the velocity spectrum 
tensor, Iiij(k, t, t'). The Galilean invariance of velocity correlators stated in equation (|4.5I) 
is most compactly expressed in Fourier-space; this is stated in the theorem below. Let 
v{K, t) be the spatial Fourier transform of v{X, r), defined by 



v{K,t) = / d^Xexp{-iK ■ X)v{X,t); [K-v{K,t)] = (5.13) 
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It is proved in the Appendix that a G-invariant Fourier-space two-point velocity corre- 
lator must be of the form 



{v^{K,T)i,*{K',T')) = {2nf6{k~k')U,j{k,t,t') 



(5.14) 



We first work out Rjm{r,t,t') and Qjmi{r,t,t') in terms of Iljm{k,t,t'). From equa- 
tion (|4.13p and (|4.14p we have 



(27r)3 (27r)3 
X exp 



d^/f d^iiT' 



2' 



(27r)3 (27r)3 
X exp 



i-iKl) {v,{K,t)v*^{K',t')) X 



i [K.X,[-,t]-K'.X, 



Substituting for Xc from equation (|4.1[) . we can write the phase 

K.X,(|,t)-K'.X,(-|,t') = (fc + fe')-^ 
where k and fc' are defined in equations (IA6[) . Therefore, 
Rj„i{r,t,t') ^ J d^kUj„,{k,t,t') exp[ik-r] 



(5.15) 



(5.16) 
(5.17) 



Qj,niir,t,t') = -i J d^A: [k - St'Snk2] Ilj,n{k,t,t') cxp[ifc.r] (5.18) 

Using equations ()4.17|) we can write the real-space integral kernels, Cjmiifjtjt') and 
Djm{r,t,t'), as 

= J d^k [nj,nik,t,t') - St'5^2'n.ji{k,t,t')]exp[ik-r] 



Cj„u {r,t,t') = Qj^i (r , i, t') + S{t - t')5m2 Qjii (r , t, t') 



= -i J d^k[ki- St'Snk2] [Ilj^{k,t,t') + S{t - t')S,n2nji{k,t,t')]exp[ik- r] 

(5.19) 
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Using equations (|5.9p we can express the Fourier-space integral kernels, Ijmi{k, t, t') and 
J jTn^k^tjt ), as 



Thus, we have expressed the integral kernels in terms of the velocity spectrum tensor, 
Tijm{k,t,t'), which is the fundamental dynamical quantity that needs to be calculated 
before the integro-differential equation for the mean magnetic field can be solved. 

6. Conclusions 

We have formulated the problem of large-scale kinematic dynamo action due to tur- 
bulence in the presence of a linear shear flow, in the limit of small magnetic Reynolds 
number (Rm) but arbitrary fluid Reynolds number. The mean-field theory we present 
is non perturbative in the shear parameter, and makes systematic use of the shearing 
coordinate transformation and the Galilean invariance of the linear shear flow. Using 
Reynolds averaging, we split the magnetic field into mean and fluctuating components. 
The mean magnetic field is driven by the Curl of the mean EMF, which in turn must be 
determined in terms of the statistics of the velocity fluctuations. In order to do this it is 
necessary to determine the magnetic fluctuations in terms of the mean magnetic field and 
the velocity fiuctuations. So we develop the equation for the fluctuations perturbatively 
in the small parameter, Rm. Using the shearing coordinate transformation, we make an 
explicit calculation of the resistive Green's function for the linear shear flow. From the 
perturbative scheme it is clear that the fluctuations can be determined to any order in 
Rm. Here we determine the magnetic fluctuations and the mean EMF to lowest order 
in Rm. The transport coefficients are given in general form in terms of the two-point 
correlators of the velocity fluctuations. At this point we make use of Galilean invari- 
ance, which is a fundamental symmetry of the problem. For Galilean invariant velocity 
statistics we prove that the transport coefficients, although space-dependent, possess the 
property of translational invariance in sheared coordinate space. An explicit expression 
for the Galilean-invariant mean EMF is derived. 

We put together all the results in § 5 by deriving the integro-differential equation 
governing the time evolution of the mean magnetic field. Some important properties of 
this equation are the following: 

(a) Velocity fluctuations contribute to two different kinds of terms, the "C" and "D" 
terms, in which first and second spatial derivatives of the mean magnetic field, respec- 
tively, appear inside the spacetime integrals. 

(6) The "C" terms are a generalization to the case of shear, of the "a" term famil- 
iar from mean-field electrodynamics in the absence of shear. However, they can also 
contribute to "magnetic diffusion" ; see discussion below. Likewise, the "D" terms are a 
generalization to the case of shear, of the "magnetic diffusion" term familiar from mean- 
field electrodynamics in the absence of shear. It must be noted that the generalization is 
non perturbative in the shear strength. 

(c) In the mean-field induction equation, the "D" terms are of such a form that: (i) 




d?k'Gr,{k-k\t,t') [nJr,^{k\t,t')^St'5r,^2^Jl{k\t,t')\ 






(5.20) 
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the equations for Bi or B3 involve only Bi or B3, respectively; (ii) the equation for B2 
involves both Bi and B2 ■ Therefore, to lowest order in Rm but to all orders in the shear 

strength, the "D" terms cannot give rise to a shear-current assist ed dynamo effect. 

(d) I n the formal limit of zero resistivity, the quasilinear theory of lSridhar fc Subramaniaii 



(|2009allbh IS recovered. In this case, the "C" terms vanish when the velocity field is non 



helical. However, this may not be the case when the resitivity is non zero. Whether the 
"C" terms give rise to such a shear~current~type effect depends on the form of the veloc- 
ity correlators, which will be strongly affected by shear and highly anisotropic; hence it is 
difficult to guess their tensorial forms a priori and it is necessary to develop a dynamical 
theory of velocity correlators - see below for further discussion. 

(e) Sheared Fourier space is the natural setting for the mean magnetic field; the normal 
modes of the theory are a set of shearing waves, labelled by their sheared wavevectors. 

(/) We prove a theorem (in the Appendix) on the form of the two-point velocity 
correlator in Fourier space; the velocity spectrum tensor and its general properties are 
discussed. The integral kernels are expressed in terms of the velocity spectrum tensor, 
which is the fundamental dynamical quantity that needs to be specified to complete the 
integro-differential equation description of the time evolution of the mean magnetic field. 

The physical meaning of the "C" and "D" terms becomes clear in the limit of a slowly 
varying magnetic f ield, when the integro -differential equation reduces to a partial dif- 
ferential equation ( Singh fc Sridhaill2010() . Then we encounter the well-known a-effect 



and turbulent magnetic diffusion (77), albeit in tensorial form. The "C" terms alone con- 
tributes to a, whereas both "C" and "D" terms contribute to magnetic diffusion. When 
the velocity field is non helical, the velocity spectrum tensor is real, and the tensorial 
a coefficient vanishes; this result is true for arbitrary values of the shear parameter. 
The "C" terms can, in principle, contribute to a shear-current effect, through the off- 
diagonal components of the diffusivity tensor (which couple the streamwise component 
of the mean magnetic field with the cross-stream components). It turns out that these 
off-diagonal components depend on the microscopic resistivity in such a manner that 
they vanis h when the microscopic resistivit y vanishes. This result is consistent with the 
results of Sridhar fc SubramanianI (|2009al lbl). To deal with the case when the micro- 



scopic resistivity does not vanish, it is necessar y to provide our kinema tic development 
with a dynamical model for the velocity field. ISingh fc Sridh^ (I2OIOI) show that, for 



forced non helical driving at low fluid Reynolds number, the sign of the off-diagonal 
terms of the diffusivity tensor does not favour the shear-current effec t. This co nclusion 
agrees with those reported in iRadler fc Stepanov (2006); Riidi ger fc Kitchatinov. (2006) : 
Brandenburg et al. ( 20081 ). even if our results are limited to low Reynolds numbers. If 



we seek a different explanation for the dynamo action seen in numerical simulations, 
the "fluctuating a-effect" still remains a promising candidate, a itself is described by 
second-order velocity correlators, so to describe fiuctuatons of a, it is necessary to deal 
with either fourth-order velocity correlators or products of two second-order velocity 
correlators. This requires extending our perturbative calculations by at least two higher 
orders, a task which, while tractable, is beyond the scope of the present investigation. 

We thank Axel Brandenburg, Karl-Heinz Radler and Matthias Rheinhardt for valuable 
comments. 
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Appendix A. A theorem on Galilean invariant velocity correlators 

Theorem 1. A G-invariant Fourier-space two-point velocity correlator must be of 
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the form 

{v,{K,t)v*{K',t')) = {2TT)H{k-k')n,j{k,t,t') (Al) 
where Ily is the velocity spectrum tensor, which must possess the following properties: 

n,,(fc,t,t') = n*^{-k,t,t') = Uj,{-k,t',t) 

K,Il,j{k,t,t') = (k,- StSak2)m,{k,t,t') = 

K'jll,j{k,t,t') = [kj - St' 5jik2)n,j{k,t,t') =0 (A 2) 

Proof. The velocity correlator in Fourier-space is 

{v^{K,t)v*{K\t')) ^ J d'^Xd'^X' e^p[i{K'-X' -K-X)]{v.,{X,t)vj{X',t')) 

= J d^X d^X' exp [i {K'- X' -K-X)] x 

X {u,{X + X,{^,t),t)v,{X' + X,{i,T'),T')) (A3) 

where equation (j4.5p has been used. Using new dummy variables of integration, X 
X — Xc{$,t) and X' — >■ X' — Xc{$,t'), we write 

{UK,r)v*{K\T')) = exp[i (K- r) - K'- r'))] x 

X J d^Xd^X'exp[i {K'-X' - K-X)] {v^{X ,t) Vj{X' ,t')) 

= exp[i {K.X,{tT)-K'.X,{tT'))] X {i,{K,T)i*{K',T')) 

(A 4) 

Comparing the left and right sides, we conclude that the phase, [K- Xc{^, t) K' ■ Xc(£, r')] , 
must vanish. Substituting for Xc from equation ()4.ip . the condition of zero phase implies 
that 

(fci - fc;)^! + (fc2 - + (fca - =0 (A 5) 

where k = (fci, ^2, fca) and fc' = (fej^, fcg) are sheared wavevectors which are related to 
K and K' through the Fourier-space shearing transformation 

fci = K1 + STK2, k2 ^ K2, ks = K3, t ^ T 

k[ = K[ + St'K'^, K'^, fcs = ^3. t' ^ ^' (A 6) 

Since equation (jA Sp must be valid for arbitrary (^1,^2,^3)1 must have A; = fe'. In 
other words, the G-invariant Fourier-space velocity correlator must be of the general 
form stated in equation (jA ip . Moreover the listed properties of the velocity spectrum 
tensor, Ily, given in equations (lA 2[) follow from the reality of v{X,t), symmetry with 
respect to simultaneous interchange of («, j) and {t,t'), and incompressibility. □ 



